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Topologically quantum error corrected logical gates are complex. Chains ol errors can form in 
space and time and diagonally in spacetime. It is highly nontrivial to determine whether a given 
logical gate is free of low weight combinations of errors leading to failure. We report a new tool 
Nestcheck capable of analyzing an arbitrary topological computation and determining the minimum 
number of errors required to cause failure. 



I. INTRODUCTION 



Constructing a quantum device for simulating quantum 
systems was first suggested by Feynman in 1982 pQ. The 
specific device he proposed, however, was not software 
programmable. The first proposal for a programmable 
quantum computer was given by Deutsch in 1985 The 
paradigm of quantum computation was thus born, but it 
had yet to acquire motivation. This came in 1994 with 
the advent of Shor's factoring algorithm [5] . 

A concern raised early within the quantum computation 
community was that of its infeasability due to decoherence 
[4]. Although this concern has yet to be eliminated, it 
has been considerably lessened. With the invention of 
quantum error correction in 1995 [5J [5] and the proof of 
the threshold theorem in 1996 [5] there is now a focus 
on constructing quantum error correction codes. 

Topological approaches to quantum error correction are 
particularly promising [9lfl3]. One such class is surface 
codes. The specific surface code we consider in this paper 
is explained well in [14] . The elementary systems of 
the surface code are physical qubits (two level quantum 
systems) |15j . and the information elements are logical 
qubits (encoded over the physical qubits). This code 
allows efficient implementation of quantum algorithms 
and its physical requirements are modest: a 2D lattice 
of qubits, nearest neighbor interactions, parallelizability, 
and gate error rates around 1% [11)1 117j . 

During computation, various gates are continually per- 
formed on the physical qubits with each gate having an 
error rate. These errors can form chains in space and 
time and diagonally in spacetime. Determining the error 
correction strength of a proposal for a surface code logical 
gate, which may consist of many of physical operations, 
can thus be somewhat complex. It can only reliably be 
accomplished via directly simulating the logical gate. 

We report a new tool for this task: Nestcheck. Using 
Nestcheck we are able to analyse an arbitrary surface 
code logical gate and determine the minimum number of 
errors required to cause failure. 



II. THE SURFACE CODE 

Group theoretic concepts provide a convenient formal- 
ism for describing both the logical qubits and gates of the 
surface code [TBI [lj|] . Homology also lends itself to our 
subject [12 [20], but encasing our discussion in the terms 
of homology might render it inaccessible. The following 
is an extremely brief account of the surface code. 

Consider a set of n qubits, Q. Decoherence of | Q) is 
due to Q interacting with an unmonitored bath, B. This 
may be written U(\Q) \B}), where U is a unitary operator. 
The effective reduced density operator is then 

Pq = tv B (U\Q) \B) (Q\ (B\ tf) = Y,E\Q) (Q\ E\ 

EES 

(1) 

where £ are Krauss operators. We call £ the errors and 
often relax the Krauss normalization condition |21j . By 
Eq. [IJ and the linearity of quantum mechanics, in order 
to correct decoherence we need only correct a basis for 
£ [22 . A convenient choice is the complete pauli basis, 
Q = {l,-l,i,-l}x{I,X,Y,Z}® n . This is a group under 
multiplication. 

The stabilizer of | Q) is the set of operators S C G such 
that S | Q) = \ Q). Every element of Q either commutes or 
anticommutes with each other element. Consider some 
S G S (a generating set for S) and applying an operator 
G G Q: If [G, S] = then SG\Q)=G\Q) and measuring 
S yields 1, whereas if {G, S} = then SG\Q) = -G\Q) 
and measuring S yields —1. 

The surface code is a stabilizer code: we enforce a cer- 
tain set of operators to be S, a set of stabilizer generators. 
We construct it such that every sufficiently low-weight 
and distinct error anticommutes with a unique subset of 
the phase factor 1 version of S. At regular intervals we 
measure the phase factor 1 version of S, each called a 
round of error correction. This yields an ordered set of Is 
and —Is, which we call a syndrome. We call an element 
of a syndrome a syndrome result. If one syndrome differs 
from another then there has been either measurement or 
physical qubit error on some number of physical qubits. 
Assuming only those low-weight errors mentioned above, 
we can determine from syndrome changes the errors that 
caused them. 



2 



IQi) 

IQ 2 ) 

\Qa) 
10} 

\Qi) 
IQ 2 ) 
IQ3) 
IQ4) 













Mx 


< 


> ( 


I 1 


I ( 


1 — 


— ( 


)— 

r 


y- 








r 








) 

( 


) 



(a) 



— ( 










M z 












1 


1 










> 

* 


1 






( 


1 



(b) 



FIG. 1: Quantum circuits for measuring (a) XX XX and (b) 
ZZZZ. 



Concerning the surface code specifically, elements of 
the phase factor 1 version of S have the form XX XX, 
XXX and XX, and ZZZZ, ZZZ and ZZ. Measuring 
the phase factor 1 version of S requires ancilla qubits , 
which we call syndrome qubits. See Fig. [T] We call the 
qubits in which logical information is stored data qubits. 
We separate S into two generating sets: primal, the X- 
operators, and dual, the Z-operators. Primal and dual 
stabilizer generators are used to correct Z- and X-errors 
respectively. This is sufficient because an /-error is trivial 
and Y = aXZ. 

Consider Fig. [2] This is a surface code plate protecting 
one logical qubit. If a Z-error acts on the center data 
qubit then the primal stabilizer generators adjacent to it 
reverse phase causing two syndrome result changes. In 
larger surface code plates chains of errors can form, the 
stabilizer generators at the ends of which reverse phase. 
Chains can undetectably connect to boundaries. In Fig. [2] 
the shortest undetectable chain of errors has weight three. 
For example, an X-error on the center-left, center and 
center-right data qubit. We say 'distance three', write 
d — 3, and this is a measure of the strength of a surface 
code logical gate. 

The normalizer of S is the set of operators AT such 
that (VAT G AO [N, S] = 0. This is the set of logical 
operators. By defining particular physical operators as 
particular logical operators we define basis logical states. 




FIG. 2: A surface code plate protecting one logical qubit. 
Filled and empty circles denote syndrome and data qubits re- 
spectively. Darkly and lightly shaded shapes represent primal 
(X) and dual (Z) stabilizer generators respectively. 



For example, in Fig. [2] a natural choice for X^ (logical 
X) is any chain of AT-operators from the left edge to the 
right that commutes with all dual (Z) stabilizers. 

The specific error correction of the surface code pro- 
ceeds via pairing syndrome results to each other and to 
boundaries. This is because every correctible error may 
be decomposed into a set of error chains each of which 
causes either two syndrome result changes, or one due to 
the error chain connecting to a boundary. Note that an 
error chain that connects to two boundaries commutes S 
and is hence uncorrectible. According that shorter error 
chains are more likely than longer ones, syndrome results 
should be paired using a minimum weight matching al- 
gorithm. Autotune, the tool we use for this task and 
which is breifly discussed in section [rET] uses Edmond's 
minimum weight perfect matching algorithm [171 124H26) . 



III. NESTCHECK 

A single round of measuring the stabilizers of Fig. [2] 
does not change the encoded data. In other words, Jl is 
applied. Less trivial logical gates involve rounds of error 
correction during which data qubits may be initialized, 
measured, swapped, and so on. Furthermore, we are not 
so fortunate that errors can only occur on data qubits 
between rounds of error correction: errors can occur on 
any qubit at any time, even measurements are prone to 
error. 

In order to determine the error correction strength of 
a surface code logical gate we first construct a primal 
and a dual graph. The vertices of such a graph are the 
spacetime locations of potential syndrome result changes 
and the edges correspond to potential connections be- 
tween syndrome result changes due to a single error. We 
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Command Figure 



Action 



ACTIVE 

INACTIVE 
ADD_X 

ADD_Z 
TRIM_X 
TRIM_Z 
HAD 

SWAP_NW 
SWAP_NE 
SWAP_SW 
SWAP_SE 
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Either SWAP-application due to a diagonally neighboring data qubit having a SWAP_NW, SWAP_NE, 
SWAP_SW or SWAP_SE command, or A(X)-application due to measurement of neighboring stabilizer 

generators. 

Consider data qubit as nonexistent. 

|+) -initialization followed by A(A)-application due to measurement of neighboring stabilizer 

generators. 

| (^-initialization followed by A(X)-application due to measurement of neighboring stabilizer 

generators. 

A(A)-application due to measurement of neighboring stabilizer generators followed by 

X-basis-measurement . 

A(X )-application due to measurement of neighboring stabilizer generators followed by 

Z-basis-measurement . 

A(A)-application due to measurement of neighboring stabilizer generators followed by //-application. 

SWAP-application of this data qubit and the syndrome qubit to the north followed by 
SWAP-application of the syndrome qubit to the north and the data qubit to the north-west. 

SWAP-application of this data qubit and the syndrome qubit to the north followed by 
SWAP-application of the syndrome qubit to the north and the data qubit to the north-east. 

SWAP-application of this data qubit and the syndrome qubit to the south followed by 
SWAP-application of the syndrome qubit to the south and the data qubit to the south-west. 

SWAP-application of this data qubit and the syndrome qubit to the south followed by 
SWAP-application of the syndrome qubit to the south and the data qubit to the south-east. 



(a) 



Command Figure Action 
ACTIVE • Either use syndrome qubit as an ancilla for stabilizer generator measurement where the data qubits 

included are those neighboring data qubits that have neither the command INACTIVE, nor SWAP_NW, 
SWAP_NE, SWAP_SW nor SWAP_SE, or SWAP-application due to a neighboring data qubit having one of 

the latter four commands. 

INACTIVE Consider syndrome qubit as nonexistent. 

(b) 

TABLE I: (a) Possible data qubit commands, (b) Possible syndrome qubit commands. Note that in figures stabilizers are 
represented as per Fig. [2] for convenience. This is not required in our set of commands because in any lattice of qubits we index 
the qubits with i and j such that if i + j mod 2 = we have a syndrome qubit and if also i mod 2 = we have a dual syndrome 
qubit. 



construct this graph using the tool Autotune [27] , 

In the Autotune lexicon, our graph of vertices and edges 
is a nest of balls and sticks. From section Ql] the data 
of error correction are the various measurements made 
during computation. Accordingly, Autotune creates a 
nest by processing many sets of measurements. A set 
is located in spacetime and has the property that if its 
measurements multiply to —1, as opposed to 1, then its 
location is that of a syndrome result change. This im- 
plies either measurement error or phase-reversal of the 
associated stabilizer generator and is called a detection 
event. In a nutshell: qubits are simulated, during which 
course all possible errors are generated and propagated, 



measurements of these qubits are placed in sets, all pos- 
sible detection events are formed to create balls, and 
pairs of detection events generated by the same error are 
connected by a stick. We thus obtain our nest of balls 
and sticks. Note that the proximity of sets to temporal 
and spatial boundaries must be specified by the user so 
that unpaired detection events generate a stick to this 
boundary specifically. 

The user must place every measurement in either two 
sets, or one set and a boundary set. Deciding when to 
create a set and in which set a particular measurement 
must be placed depends on the physical gate sequence 
the user is applying. This physical gate sequence can 
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Figure 



Description 



/L-application. 
1 0)-initialization. 
^-application, 
^-basis-measurement . 
A(X)-application (left qubit is the control). 

SWAP-application. 



t 
I 
\ 

H 



l_ _' A set, whose measurements are those enclosed. 

'*...<• The measurement enclosed is in a boundary set. 

TABLE II: Notation scheme for parts of figures [3] to [7] that 
display physical gate sequences specified by frames. 



be rather complex: a surface code computation typically 
involves introduction, deformation and elimination of 
boundaries, measurements in various bases, and single 
and double physical qubit gates. We desire a short list of 
commands to be given locally to physical qubits such that 
an arbitrary surface code computation can be specified. 
Table UJ details the list we use. Note that initializations 
specific to state-injection have not been included, but will 
be in future work. 

We construct a surface code computation by writing 
a program. A program is a series of frames. A frame 
specifies a single round of error correction and is a list of 
physical qubits each with two data: a command and a 
boundary specification. We require boundary specifica- 
tions for while Autotune can discern when a boundary set 
is required, it cannot discern which boundary set. Figs. [3] 
to [7] display small programs and elucidate how these spec- 
ify gate sequences and sets. From these figures the logic 
of creating and placing measurements in sets can in part 
be infered. We expand on this below. 

There is some arbitrariness in how we use frames to 
specify sets. In sets containing two measurements, the 
measurements often have different times. The two mea- 
surements are specified by different frames and so the 
set is specified by these two frames conjointly. It is thus 
somewhat arbitrary which frame we choose to specify the 
boundary set that this set may connect to. We choose 
the frame executed at a later time. 

It is convenient to introduce some terminology for dis- 
cussing frames: If every data qubit has either the com- 
mand INACTIVE or ACTIVE we call the frame an I L -frame. 
Otherwise, we call the frame by the commands of the data 
qubits that are not INACTIVE or ACTIVE. For example, if 
some data qubit has the command HAD we call the frame 
a HAD-frame. Obviously, this is not the best terminology 






FIG. 3: A program consisting of two /L-frames. (a) Represen- 
tation of the program frames: Read upwards. See table [T] for 
the notation scheme, (b) Isometric view with upwards world 
lines: See table ITT] for the notation scheme. 



for frames in general: what do we call a frame in which 
some qubits have the command ADD_X while others have 
the command ADD_Z? Nevertheless, it serves us well in 
the discussion of boundary sets below. 

In order to handle boundary sets, nonrepeating Jl- 
frames must often intersperse the frames of a program. 
By 'nonrepeating /L-frame' we mean an JL-frame in which 
more boundaries are specified than necessary were we 
to continually execute it. This is best illustrated by 
examining programs in which various frames are followed 
by /L-frames. There are four cases: 

(i) An JL-frame followed by I^-ir ames: Obviously, of 
the JL-frames following the IL-frame none needs to be 
nonrepeating. See Fig. [3] 

(ii) An ADD_X- or ADD_Z-frame followed by 7 L -frames: 
Consider specifically the ADD_X- frame shown in Fig. |4|d). 
Syndrome qubits are activated in this frame. Since we are 
|+)-initialising data qubits, the first syndrome result of a 
dual syndrome qubit activated in this frame is random 
and must be placed in a boundary set and a set that 
connects to this boundary. The user must specify which 
boundary set in the ADD_X-frame, and while Autotune 
can discern that the connection mentioned is required it 
cannot discern which boundary set must be connected 
to. The user must therefore specify the connection in the 
directly following I^-iiame. Such boundary information is 
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1+) 1+) 




'A a 

(c) (d) 

FIG. 4: Each pair (a)-(c) and (b)-(d) is a program consisting 
of a an iL-frame followed by a ADD_X-frame. (a) and (b) are 
isometric views with upwards world lines: See table [IT] for 
the notation scheme, (c) and (d) are representations of the 
program frames: Read upwards. See table |TJ for the notation 
scheme. 



not required in a repeating iL-frame. Thus the iL-frame 
that directly follows the ADD_X- or ADD_Z-frame needs 
to be nonrepeating, all other following I^-hames can be 
repeating. 

(iii) A TRIM_X- or TRIM_Z-frame followed by J L -frames: 
Consider specifically the TRIM_X- frame in Fig.[5jd). In the 
directly following frame, syndrome qubits are deactivated. 
We cannot use the X-basis- measurements of data qubits to 
construct dual (Z) syndrome results. The final syndrome 
result of a dual syndrome qubit deactivated in the directly 
following IL-frame must hence be placed in a boundary 
set and a set that connects to this boundary. Both of 
these sets are specified by the TRIM_X-frame. Thus none 



i o i 




FIG. 5: Each pair (a)-(c) and (b)-(d) is a program consisting 
of an TRIM_X-frame followed by a 7L-frame. (a) and (b) are 
isometric views with upwards world lines: See table [TT] for 
the notation scheme, (c) and (d) are representations of the 
program frames: Read upwards. See table [TJ for the notation 
scheme. 



of the J L -frames that follows the TRIM_X- or TRIM_Z-frame 
needs to be nonrepeating. 

(iv) A HAD- frame followed by a SWAP_NW-, SWAP_NE-, 
SWAP_SW- or SWAP_SE-frame followed by 7 L -frames: Con- 
sider specifically the SWAP_NW- frame in Fig. [6^b). In 
the directly following frame, syndrome qubits are acti- 
vated. Ideally, the first syndrome result of each syndrome 
qubit activated in the directly following frame should be 
matched with that of the previous HAD-frame diagonally 
offset from it. This would involve, however, converting a 
primal nest into a dual nest and vice versa, which would 
in practive be quite a bit of effort code. Perhaps in future 
research we shall concern ourselves with this, but for the 
moment we use a shortcut. We simply consider activated 
syndromes as entirely new. The first syndrome result of 
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FIG. 6: A program consisting of a HAD-frame followed by 
a SWAP_NW-frame followed by an /L-frame. Note that these 
frames are to be understood as possibly extending an arbitrary 
distance in the north, west, east and south directions, lest data 



FIG. 7: A program consisting of a HAD-frame followed by 
a SWAP_NW-frame followed by an /L-frame. Note that these 
frames are to be understood as possibly extending an arbitrary 
distance in the north, west, east and south directions, lest data 



qubits appear from and dissapear to nowhere (see Fig. 11 d) 
of appendix |E] for an example of a complete swapping), (a) 
Representation of the program frames: Read upwards. See 
table [I] for the notation scheme, (b) Isometric view with 
upwards world lines: See table III] for the notation scheme. 



qubits appear from and dissapear to nowhere (see Fig. 11 d) 
of appendix [E] for an example of a complete swapping), (a) 
Representation of the program frames: Read upwards. See 
table [I] for the notation scheme, (b) Isometric view with 
upwards world lines: See table III] for the notation scheme. 



each syndrome qubit activated in the directly following 
frame is then considered random and must be placed in a 
boundary set and a set that connects to this boundary. 
This situation is similar to case (i) but one frame on- 
wards. Thus two of the iL-frames that directly follow the 
SWAP_NW-, SWAP_NE-, SWAP_SW- or SWAP_SE-frame need to 
be nonrepeating. 

A further note is necessary on case (iv) : The HAD-frame 
requires more boundaries than one might naively expect to 



be specified. Consider specifically the SWAP_NW-frame in 
Fig. [7](b). In this frame syndrome qubits are deactivated. 
Since we are using a shortcut, the final syndrome result 
of each syndrome qubit deactivated in this frame must 
be placed in a boundary set and a set that connects to 
this set. Both of these sets are specified by the previous 
frame, the HAD-frame. Thus the HAD-frame must include 
boundary information to this effect. 

An indispensible tool used in debugging the logic for 
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FIG. 8: Color online. Blender model of a section of the 
primal nest for d = 3 -Hl as presented in appendix [B] Thicker 
cylinders (sticks) represent more probable connections between 
potential syndrome result changes (balls) due to a single error. 
Red cylinders represent sticks between sets, differently colored 
cylinders represent sticks between a set and a boundary set. 
In this nest section we see trimming and the beginning of 
swapping. There are more sets than usual in a trimming 
time-slice, hence the offset cylinder intersection points. 



handling sets is a simple Blender visualizer we have devel- 
oped [57] . This tool creates a 3D blender model of a nest 
so that we can see whether the correct boundaries have 
been specified and the correct connections exist. Fig. [8] is 
such a model of a section of the primal nest for d = 3 
as set out in appendix [Bj 

With the primal and dual nests built for a particular 
surface code logical gate we then need to find the short- 
est chain between any two boundaries and the smallest 
topologically nontrivial ring. The length of the smaller of 
these two is the distance of this surface code logical gate 
and so the error correction strength. 

Finding the shortest chain between two boundaries is 
simple. We loop over pairs of boundaries and perform a 
breadth first search between them. 

As to finding the smallest topologically nontrivial ring 
the situation is not quite so undemanding and we have not 
yet implemented such a search. This is not to invalidate 
Nestcheck, however, for rings are generally easy to see. 
Nevertheless, without this ability Nestcheck cannot be 
considered complete. 



In future work we plan to search for a set of smallest 
topologically distinct and nontrivial rings. An outline of 
our prospective method is thus: The 'hole' around which 
a topologically nontrivial ring in a nest of one type can 
exist is itself a ring, or at least a closed structure [35], of 
the other type. This latter ring can be specified by the 
user. The task is then reduced to, for each input ring of 
a certain type, looping over an appropriately constructed 
set of rings in the nest of the other type and checking 
whether each interlocks with the input ring. The concep- 
tually principal part of this algorithm is a method for 
determining whether a given primal ring interlocks with 
a given dual ring. 

Nestcheck has been tested extensively on a wide variety 
of physical gate sequences. Furthermore, as an example, 
we have run Nestcheck on a proposal for d = 3 (see 
appendix EJ. Nestcheck found that this proposal is d = 3 
as purported. 



IV. CONCLUSION 

We have reported a new tool Nestcheck capable of 
analysing an arbitrary surface code computation and de- 
termining the minimum number of errors required to 
cause failure. 3D topological cluster states may also be 
analysed since the data that Nestcheck analyses is created 
using Autotune. In future research, surface code compu- 
tations consisting of many qubits shall be constructed, 
for which it is highly nontrivial to check the strength of 
error correction. Nestcheck shall therefore be of utility 
in this endeavour. Furthermore, using the framework of 
Nestcheck (programs, frames, and so on), we are able to 
easily specify a complex surface code computation thus 
enabling simulation of such to be relatively quickly and 
painlessly written. 
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been included. The final frame, Fig.[9]jc), is repeating in 
the sense of section Mil 

Appendix B: Distance three logical hadamard 



A scheme for d = 7 Hl is presented in [33] . We adapt 
this to d — 3. Besides nonrepeating I^-irames (see sec- 
tion [Tnl Figs. [TU] to [12] present this adaptation. 
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FIG. 9: Full program for simulating the surface code plate of 
Fig. [2] Read (a), (b) then (c). See table [TJ for the notation 
scheme. Arrows represent boundary specifications: solid and 
dashed for primal and dual respectively, and bent and straight 
for temporal and spatial respectively. We see three primal 
boundaries: a temporal one, one along the top and one along 
the bottom. We see two dual boundaries: one along the left and 
another along the right. Note that in general spatiotemporal 
boundary specifications are possible. 



Appendix A: Full program for a distance three plate 

Fig. [9] details an entire program for simulating the 
surface code plate of Fig. [2] Boundary specifications have 



Appendix C: Checking the logic of specific surface 
code logical gates 



There would be little point in checking the error cor- 
rection strength of a proposal for surface code logical 
gate without checking whether it is indeed the logical 
gate purported. Future research shall potentially yield 
a general tool for this task [55], but for the moment we 
resort to specific methods. 

Excepting initializations specific to state- injection, the 
physical gates of an arbitary surface code computation are 
in the Clifford group. Algorithms for efficiently simulat- 
ing Clifford group quantum circuits have been developed 
and implemented, and are available for download. Wc 
currently use Graphsim [3D]. Constructing a Graphsim 
simulation is simple: we declare a qubit register and then 
call our gate sequence. 

To view the state of the qubit register at a point during 
simulation there are two methods. We may either print 
a tabic of stabilizers, or simply create a duplicate of the 
qubit register and measure all the qubits in an appropriate 
basis. Since we are checking rather than investigating 
logic we choose the latter method. 

Take J?l as an example. Here we must check whether 
the gate maps |0), |+) and |-) to |+), |-), |0) and 
|1) respectively. Consider the case |0) H> |+): We ini- 
tialize the logical qubit to |0) and then apply the gate 
sequence of Hi, , during which course byproduct-operators 
[13] accrue. Then we measure in the X-basis each qubit 
to which Xi is nontrivial and check whether the multiple 
of these measurements with those associated with the 
appropriate byproduct-operators is 1. Success is declared 
if so, otherwise failure. 

This method of checking logic is stochastic. In order to 
reliably check using this method, then, we must apply it 
a large number of times. If failure is declared even once 
then the proposal fails entirely. 



[1] R. P. Feynman, Int. J. Theor. Phys. 21, 467 (1982). 
[2] D. Deutsch, Proc. R. Soc. Lond. A 400, 97 (1985). 



[3] P. W. Shor, in Proc. 35th Annual Symposium on Founda- 
tions of Computer Science (IEEE Computer Society Press, 



9 



Los Alamitos, CA, 1994), pp. 124-134, quant-ph/9508027. 

[4] A. Berthiaume, D. Deutsch, and R. Jozsa, in Physics and 
Computation, 1994- PhysComp'94, Proceedings, Work- 
shop on (IEEE, 1994), pp. 60-62. 

[5] P. W. Shor, Phys. Rev. A 52, R2493 (1995). 

[6] A. M. Steane, Phys. Rev. Lett. 77, 793 (1996). 

[7] E. Knill, R. Laflamme, and W. Zurek, quant-ph/9610011 

(1996) . 

[8] D. Aharonov and M. Ben-Or, Proc. ACM STOC 29, 176 

(1997) , quant-ph/9611025. 

[9] S. B. Bravyi and A. Y. Kitaev, quant-ph/9811052 (1998). 
[10] E. Dennis, A. Y. Kitaev, A. Landahl, and J. Preskill, J. 

Math. Phys. 43, 4452 (2002), quant-ph/01 10143. 
[11] R. Raussendorf, J. Harrington, and K. Goyal, Ann. Phys. 

321, 2242 (2006), quant-ph/0510135. 
[12] R. Raussendorf, J. Harrington, and K. Goyal, New J. 

Phys. 9, 199 (2007), quant-ph/0703143. 
[13] A. G. Fowler and K. Goyal, Quant. Info. Comput. 9, 721 

(2009), arXiv:0805.3202. 
[14] A. G. Fowler, M. Mariantoni, J. M. Martinis, and 

A. N. Cleland, Phys. Rev. A 86, 032324 (2012), 

arXiv:1208.0928. 
[15] B. Schumacher, Physical Review A 51, 2738 (1995). 
[16] D. S. Wang, A. G. Fowler, and L. C. L. Hollenberg, Phys. 

Rev. A 83, 020302(R) (2011), arXiv: 1009.3686. 
[17] A. G. Fowler, A. C. Whiteside, and L. C. L. Hollenberg, 

arXiv:1202.5602 (2012). 



[18] D. Gottesman, Ph.D. thesis, Caltech (1997), quant- 
ph/9705052. 

[19] D. Gottesman, Phys. Rev. A 57, 127 (1998), quant- 
ph/9702029. 

[20] R. Raussendorf and J. Harrington, Phys. Rev. Lett. 98, 

190504 (2007), quant-ph/0610082. 
[21] P. Kaye, R. Laflamme, and M. Mosca, An Introduc- 
tion to Quantum Computing (Oxford University Press, 

Great Clarendon Street, Oxford OX2 6DP, 2007), ISBN 

9780198570493. 
[22] M. A. Nielsen and I. L. Chuang, Quantum Computation 

and Quantum Information (Cambridge University Press, 

Cambridge, 2000). 
[23] D. P. DiVincenzo and P. W. Shor, Phys. Rev. Lett. 77, 

3260 (1997), quant-ph/9605031. 
[24] J. Edmonds, Canad. J. Math. 17, 449 (1965). 
[25] J. Edmonds, J. Res. Nat. Bur. Standards 69B, 125 (1965). 
[26] A. G. Fowler, A. C. Whiteside, and L. C. L. Hollenberg, 

Phys. Rev. Lett. 108, 180501 (2012), arXiv:1110.5133. 
[27] A. G. Fowler, A. C. Whiteside, A. L. Mclnnes, and A. Rab- 

bani, arXiv:1202.6111 (2012). 
[28] A. G. Fowler and S. J. Devitt, arXiv:1209.0510 (2012). 
[29] A. G. Fowler, arXiv: 1202.2639 (2012). 
[30] S. Anders and H.-J. Briegel, Phys. Rev. A 73, 022334 

(2006), quant-ph/0504117. 



10 




o.o.o o . 
.0.00.0.0.0.0 

IT TUT 




M z 



M 



M z 



M, 



W \J w 



M 



^^^^^ ^^^^^ 
O • O • 



o . T 



M z M x M x 
.0.0 o . 

0.0.0.0 

.0.0.0. 

Ml JL. 



M, 



M, 



M z W A 
0.0.0.0 

^A/ z Afx A/ x 
.0.0 o . 




• <> . <> ^o^ . ^o^ . ^o^ 



H 



H 



H 



H 



H 



H 



H 



H 



H 



H 



H 



H 



1 



\ \ \ 

. o_ . o_ • o 

\ \ 

\ \ \ 

. o_ . o_ • o 

\ \ 

\ \ \ 





■ 

o . 

o 




::: 



: : 

o . 



0.0.0.0.0 



o o 



(c) 



(d) 



FIG. 10: Frames for simulating d = 3 Hi. See table [I] for the notation scheme. 
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FIG. 11: Frames for simulating d — 3 Hi. See table [I] for the notation scheme. 
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FIG. 12: Frames for simulating d — 3 H^. See table [I] for the notation scheme. 



